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ABSTRACT 

The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be 
fully anisotropic, i.e. one has to distinguish three diffusion axes in a local, field-aligned frame. We 
reexamine the transformation for the diffusion tensor from this local to a global frame, in which the 
Parker transport equation for energetic particles is usually formulated and solved. Particularly, we 
generalize the transformation formulas to allow for an explicit choice of two principal local perpendic- 
ular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case 
of isotropic perpendicular diffusion. For the local frame, we motivate the choice of the Frenet-Serret 
trihedron which is related to the intrinsic magnetic field geometry. We directly compare the old and 
the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk 
and the Parker field. Subsequently, we examine the significance of the different formulations for the 
diffusion tensor in a standard 3D model for the modulation of galactic protons. For this we utilize 
a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker 
transport equation and present the resulting modulated spectra. The computed differential fluxes 
based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid 
for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending 
on heliocentric distance. 

Subject headings: cosmic rays — diffusion — Sun: heliosphere 



1. INTRODUCTION AND MOTIVATION 

The most important transport process for energetic 
charged particles in the heliosphere is their spatial diffu- 
sion as a consequence of their interaction with the tur- 
bulent heliospheric magnetic field (HMF). While rarely 
used in models of galactic transport, the con cept of 
anisotropic diffusion is well established (se e, e.g., Burger 
et al.|200Q| [Schlickeiser|2002l |Shalchi|2009| ) for models of 
heliospheric cosmic ray (CR) modulation. In most stud- 
ies the anisotropy refers to a difference in the diffusion 
coefficient parallel (Ky) and perpendicular (k±_) to the 
magnetic field, and the perpendicular diffusion is treated 
as being isotropic. 

Although the notion of fully anisotr o pic d iffusion is not 
new - see an early study by Jokipii (1973) considering 
for the first time anisotropic perpendicular diffusion, i.e. 
K_L1 7^ K -L2 _ it was not before the measurements made 
with the Ulysses spacecraft that this concept had to be 
used to expla in the high-latitude observations of cosmi c 
ray s, see, e.g., Jokipii et al. ( 1995 ) , Potgieter et al. ( 1997[ ) , 
and Ferreira et al.||2l)I)T| ). These studies remained largely 
phenomenological and did not attempt a rigorous inves- 
tigation of anisotropic perpendicular diffusion. 

More recently, in the context of studies of the trans- 
port of solar energe t ic par ticles in the heliospheric P arker 



Kelly et al.| f2012| have 
le diffusion tensor from 



field, |Tautz et al.] p)llj ) and 
determined the elements of t 
test particle simulations in a local, field-aligned frame. 
While the former authors find no conclusive result, the 
latter authors clearly demonstrated that the scattering 
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in the inhomogeneous Parker field can indeed induce 
anisotropic perpendicular diffusion. 

Given this phenomenological and simulation-based ev- 
idence it is important to determine the principal direc- 
tions of perpendicular diffusion in the field-aligned local 
frame, because the transformation of the diffusion ten- 
sor from a correspondingly oriented local into a global 
coordinate system determines the exact form of the ten- 
sor elements in the latter, in which the transport equa- 
tion is usually solved. This is of particular importance 
in the case of sy mmetry- fre e magnetic fields, like the so- 
called Fisk field (|Fisk|1996|). The latter is - although in a 



weaker manner than orig inally suspected (Lionello et al. 



2006 Sternal et al.||20lT ) - still a valid generalisation of 



the Parker field and takes into account a non-vanishing 
latitudinal field component. 

While it has been recognized that the use of the Fisk 
field in models of the heliospheric modulatio n of CRs re- 
quire s a re-derivation of the diffusion tensor ( |Kobylinski| 
2001||Alania|2002j [Burger et al.|2008| the formulas given 
in these papers differ from each other and are either valid 
only for the case of isotropic perpendicular diffusion (the 
former two papers) or for a specific orientation of the lo- 
cal coordinate system (the latter paper). Consequently, 
there are two open issues, namely (i) to determine which 
of these formulas are correct (see also Appendix A) and 
(ii) to generalize these results to the case of anisotropic 
perpendicular diffusion. With the present paper we ad- 
dress both issues by deriving general formulas for the 
transformation of a fully anisotropic diffusion tensor. In 
addition to establishing the appropriate description, we 
apply the new generalized formulas to a standard modu- 
lation problem in order to demonstrate the physical sig- 
nificance of the approach. 
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Figure 1. The undisturbed (right) heliospheric magnetic field 
(projected into the equato r ial (to p) and a meridional plane (bot- 
tom)) according tolParkerl d 19581) and its s tructure when field line 
random walk is included (left), taken from |Jokipii| | |2001[ l . 



2. GENERAL CONSIDERATIONS 

Anisotropic perpendicular transport can, in principle, 
result (i) from an inhomogeneous (asymmetric) mag- 
netic background field or (ii) from turbulence that is 
intrinsically non-axisymmetric with respect to t he (ho- 
mogeneous) local magnetic field direction (e.g., Wein- 
horst et al. 20081. While the latter case has been dis- 



cussed inTTie con text of energetic particle transport ( Ruf- 
folo et al. 2008) partly motivated by the observed ra- 
tios of the power in the microscale magnetic field fluc- 
tuations parallel and perpendicular to the background 
field: 5B\ X : 5B\ 2 : 5B\ = 5 : 4 : 1 (where SB ±1 is 
aligned to the latitudinal unit vector an d the normalized 
SB ±2 completes the local trihedron, see Belcher & Davis 
1971 Horbury et al.|[l995 ), recent analyses indicate that 
the p erpendicular fluctuations are probably axisymmet- 



ric dTurner et al. 2011 Wicks et al. 2012). Therefore 



we consider the first case of an inhomogeneous mag- 
netic background field to be more likely to cause fully 
anisotropic diffusion. 

If the random walk of field lines due to turbulence is 
significantly contributing to the perpendicular particle 
transport, one generally has to expect the latter to be 
anisotropic. This can be illustrated already for the sim- 
ple case that the HMF is represented by the Parker spi- 
ral, see Fig. [T] Due to the field geometry the field line 
wandering is not isotropic, neither in radial direction nor 
in heliographic latitude, resulti ng in a field line d iffusion 



coefficient depending on both ( Webb et ah] 2009 ) . 

As soon as anisotropic perpendicular diffusion occurs, 
it is necessary to determine the principal axes of the dif- 
fusion tensor in a local field-aligned frame (kl), because 
their orientation determines the tensor elements in the 
global frame (k) after a corresponding transformation 
given by 



' K_L1 K A 
— KA ^_L2 
K\\ 



(2) 



where, in general, ka denotes the drift coefficient, in- 
duced by a non-axisymmetric turbulence and by inho- 
mogeneous magnetic fields. The latter drifts can always 
be de scribed by a drift velocity vg in the transpo rt equa- 
tion ( jTautz fc Shalchi||2012| |Burger et aL|2008[ ) and are 
therefore not considered in the following. In Eq. ( 1 ) , 
analogous to the Euler angle transformation known from 
classical mechanics, the matrix A = R 3 R 2 Ri describes 
three consecutive rotations Ri with A^ 1 = A T . These 
rotations are defined by the relative orientation of the 
local and the global coordinate system with respect to 
each other. 

Due to the latitudinal structuring of the solar wind 
and, in turn, of the Parker spiral having a vanishing B§- 
component, one may argue that in that case the latitu- 
dinal direction remains a preferred one so that the local 
coordinate system could always be defined by the unit 
vectors t (along the field), e# (from a spherical polar co- 
ordinate system) and e$ xt. This, however, can obviously 
not be the c ase for symmetry-free fields like the Fisk field 
flFisk||1996| ). 

In general, the local trihedron will consist of a unit 
vector t tangential to the magnetic field and two orthog- 
onal ones, u and v, defining the remaining principal axes. 
With this notation the transformation (1) reads, for an 
arbitrary choice of this local trihedron: 



K ll = K _L1 u \ + K _L2 v i + K || t\ 

«12 = K_L1 Ui U 2 + K_L2 V\ V 2 + K|| t\ t 2 

+ K A (Ul V 2 ~ U 2 Vi) 
«13 = «-Ll Ul U 3 + Kj_ 2 Vi V 3 + K\\ h t 3 

+ K A (Ul V 3 - U 3 Vl) 
K 2 1 = Ui U 2 + Kj_ 2 Vl V 2 + K\\ t 1 t 2 

~K A (Ul V 2 - U 2 Vl) 
K 22 = K_L1 ul + Kj_ 2 v\ + K|| t\ 
K23 = U 2 U 3 + K_L2 V 2 V 3 + K\\ t 2 t 3 

+ K A (U 2 V 3 - U 3 V 2 ) 
K31 = Ul U 3 + Kj_ 2 Vl V 3 + K|| tl t 3 

-KA (Ui V 3 - U 3 Vi) 
K 3 2 = U 2 U 3 + K ±2 V 2 V 3 + K\\ t 2 t 3 

~ k a (u 2 V 3 - U 3 V 2 ) 
K33 = K±l uj + Kj_ 2 vl + K|| tj 



(3) 
(4) 

(5) 

(6) 

(7) 
(8) 

(9) 
(10) 

(11) 



where the components of t, u and v are determined in 
the global coordinate system. Consequently, the task is 
to determine the unit vectors t, u, and v for an arbitrary, 
symmetry-free magnetic field. 
Given that t he perpendicular fluctuations are proba bly 



axisymmetric (Turner et al. 2011 Wicks et al. 



2012) as 



AklA 1 



(1) 



discussed above, we assume ka = in the following. 

With this explicit formulation of the tensor elements 
we can already address issue (i) defined in Section 1: For 
the case that the perpendicular diffusion is isotropic, i.e. 
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k±i = k±2, the formulas given by Burger et al. 
see Appendix A, are identical to Eg. ( 3) to (11) 



42008] ), 
so that 



thei r correction of the results found by Kobylinski (2001 1 
and Alania ( 2002 ) and, in turn, their subsequent analysis 
are validated! We emphasize, however, that neither of 
these formulations (involving only two rotation angles) 
allow to define explicitly the perpendicular diffusion axes, 
which are nessesary to treat anisotropic diffusion in the 
most general form. 

3. THE CHOICE OF THE LOCAL COORDINATE SYSTEM 

In the absence of symmetries there remain two distin- 
guished local directions that, at a given location within 
an arbitrary magnetic field, are related to its curvature k 
and torsion r and are called the normal and the binormal 
direction. They can be defined with the corresponding 
normal and binormal unit vectors, respectively. Together 
with the tangential unit vector, they constitute a lo- 
cal orthogonal trihedron fulfi lling the (k- and r-d efining) 
Frenet-Serret relations (e.g., Marris fc Passman||1969 ): 

(12) 
(13) 
(14) 



(f- V)t = kn 
(t ■ V)ft=-kt + rb 
(f- V)b=-rn 



If no other diffusion axes are preferred by any process, the 
Frenet-Serret System constituted by the above definition 
of t, ft, and b is the most natural choice, i.e. u = n and 

v = b in Eqs. (3) to (11). 

The transformation of the local diffusion tensor into 
a global coordinate system according to these equations 
thus requires knowledge of the dependence of the Frenet- 
Serret vectors on a given (non-homogeneous) magnetic 
field B. Evidently, the required relations are 

t = B/\B\ (15) 

n=(t-V)t/k (16) 

b = txn (17) 

This trihedron can, of course, only be established for a 
spatially non-homogeneous field, but this (weak) condi- 
tion is fulfilled in most cases of interest. If there would 
exist a region where the field is homogeneous, the choice 
of the vectors n and b is arbitrary (signifying isotropic 
perpendicular diffusion) unless no other preferential di- 
rections unrelated to the field geometry can be specified. 
Other principal directions unrelated to the large-scale ge- 
ometry of the field could, for example, arise from non- 
axisymmetric turbulence. The above formulas (3) to (11) 
remain unaffected, however: One only needs to specify 

the appropriate vectors t, ft and b. 

In the following we illustrate the procedure with the ex- 
ample of the well-studied HMF. We quantitatively com- 
pare the new tensor with the 'traditional' one, which 
is only valid for isotropic perpendicular diffusion. This 
comparison reveals that a study of fully anisotropic tur- 
bulent diffusion within more complicated fields - like 



the much-discussed heliospheric Fisk field (Burger et al. 
2008)|Sternal et al.|2011||Fisk|1996[[B"urKer OWe]iW04j ) 



or co mplex galactic magnetic fields (Kuzmaikin et al 
1988 Beck et aL||1996 ) - has to be pertormed with even 
more caution than thought before. 




Figure 2. The hybrid Fisk and the Parker field illustrated by red 
and black field lines, respectively. The two local trihedrons for the 
Parker field are indicated with the orange and blue (Frenet-Serret) 
as well as the yellow and light blue (traditional) lines. Note that 
the traditional trihedron is always aligned to the Parker spiral cone 
of constant 1? while for the Frenet-Serret trihedron one axis (the 
K ±2-binormal axis, orange) is nearly parallel to the z-direction. In 
the ecliptic both coordinate systems coincide by definition. All 
distances arc in units of AU. 



4. AN EXAMPLE FOR THE NEW DIFFUSION TENSOR 

4.1. The heliospheric magnetic field 

An analytical representation of the HMF, which is 
referred to as the h ybrid Fisk field, can be found in 
Sternal et al. (2011). For a constant solar wind speed 
400 km/ s J the HMF is represented, using spheri- 



cal polar coordinates, by the following formulation: 



B r = A B e 



Bfl — B r 



B u 



dr) 



T 



* • /I* ■ * 

— ui sin p sin if 
suvd{oj* cos f3* 



Or. 



(18) 
(19) 
(20) 



(oj* sin j3* sin ■&) cos ip* } 



with 



where 



FsW = 
' [tanh((5 p #) + tanh(5 p (# • 
-tanh(5 e (tf-tf 6 ))] 2 


[tanh(5 p i?) + tanh(5 p ($ - 
-tanh((5 e (i?-7r + t? b ))] ; 



(21) 



7T)) 



7T)) 



< < A 

fib < < 7T - ti b 
7T — < d < 7T 



is the transition function introduced by |Burger et al.| 
(2008). In the case F s = the HMF reduces to the 
standard Parker spiral magnetic field. For a quantitative 
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comparison of different HMF configurations, see |Scherer| 

in Eq. (18), B e denotes the magnetic field strength 
at r e = 1 AU, r is the solar radius, and f2© = 2.9 • 
10 -6 Hz is the averaged solar rotation frequency. The 



constant A = ±1 in Eq. (18) indicates the different field 



directions in the northern and southern hemisphere. The 
values for the angle between the rotational and the so- 
called virtual axes of the Sun j3 = 12° and the differential 
rotatio n rate ui = fi© /4 are taken from Sternal et al. 
(2011). The parameters <L = 5 and S c = 5 determine 
the respective contributions of the Fisk- and Parker field 
above the poles and in the ecliptic while = 80° is 
the cutoff colatitude for the Fisk-field influence. In the 
following, we consider two cases, a pure Parker field (i.e. 
setting F s = in Eqs. (19) and (20)) and the hybrid Fisk 
field with F s from Eq. (21). Both fields are illustrated 
by exemplary field lines in Fig. [2] 

4.2. The local diffusion tensor elements 
The elements of the local diffusion tensor are chosen 



following the approach in Reinecke et al. (1993), i.e. as 



=*'<•<' l;UUJ (22) 

— 4£Ktr <23 » 

K±2 = foil (24) 

where (3 — v/c is the particle speed normalized to the 
speed of light, p is the particle momentum with the nor- 
malisation constant p = 1 GeV/c and B is the magni- 
tude of the magnetic field. The scaling exponents have 
the values an = 0.75 and a± = 0.97. The parallel diffu- 
sion constant is «m = 0.9- 10 22 cm 2 /s while k± = 0.1kii . 
The anisotropy in perpendicular diffusion is assumed to 
be soley determined by the factor £ which is set to 2 for 
the following discussion. This is still a moderate choice 



compared to the findings of, e.g., Potgieter et al. (1997). 

Although these empirical formulas for the locardiffu- 
sion coefficients are not directly related to the turbulence 
evolution in the heliosphere and more sophisticated the- 
oretical models for the corresponding mean free paths in 
parallel and perpendicular direction exist, they are still a 
good approximation as can be seen in the following. The 
result from q uasilinea r theo ry (QLT) for the mean free 
path (see e.g. Shalchi (2009)) is given by 



A 



(QLT) 



3Z ; 



slab 



16nC(u) \5B. 



B 



'slab 



R 



2-2v 



{l-v)(2-v) 



R 



2// 



with 



C(u) = 



1 I>) 

2y/nT(v- 1/2) 



(25) 



(26) 



where T(x) is the gamma function, 2v — 5/3 is the 
inertial range spectral index, R — Rh/hiab is the di- 
mensionless rigidity, and Rl = pc/(\q\B) is the particle 
Larmor radius. If one scales the bendover scale of slab 
turbulence as l s i a b = 0.03 p 0,5 (where p is the heliocen- 
tric distance in AU) and the slab turbulence variance 



aS dB slab 



B't 



-2.15 



the radial dependence of the lo- 



cal tensor ele men ts m atc hes well with the approxima- 
tive formulas (22) to (24) as shown in Fig. [3 It is in- 
teresting to note that these scalings are sh fnlar to the 
assumptions made in |Burger et al.| (2008). They use 
the same radial dependence tor l s i a b (their exponent of 



lab 



m in = 32 p U 5 is a typing error, private com- 
munication with the authors) and an exponent of —2.5 
for the slab turbulence variance o~B^ lab , which is slightly 
larger. Similar arguments can be made for the perpen- 
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Figure 3. The dependence of the local and global tensor elements 
on heliocentric distance in the eclipti c plane f or the Parker field. 
The local elements from the formulas l |22| -( [24| ) are shown as solid 
red (k||), green (rej_x), an d blue (kx2~) curves. The results for 

the parallel diffusion coefficent kJ | QLT) = 1 /3»;Aj | QLT) (EqEiJ and 

fNLGCI / I L 

the perpendicular coefficient from re^ ' (Eq. 27| are drawn as 
dashed red and green lines, respectively, while the dashed blue 
curve is just scaled as ^k^^ GC ' with § = 2. The overlayed, color- 
matched symbolds show the nearly perfect alignment of the K rr 
(green, x), K$fi (blue, •) and k vv (red, +) global diagonal tensor 
elements in the ecliptic, due to the Parker field structure. All other 
tensor elements are almost indistinguishable from zero, as indicated 
by the remaining symbols. 

dicular diffusion. Employing for the perpendicular diffu- 
sion coefficient the result of the nonlinear g uidin g center 
(NLGC) the ory of 
et al.l d2004b, name 



Matthaeus et al. (2003) and Shalchi 



(NLGC) _ 



a v 



y 



r(i//2 + i) 



2V3V'T>/2 + l/2)" 2Z5 ~ 5 



SBln 



-, 2/3 



.1/3 



(27) 



(see formula (15) in Burger et al. (2008)) with the con- 
stant a = l/v3- Scaling again the 2D turbulence corre- 
lation length I2D with p 5 and the turbulence variance 

yields resu lts simi- 



(5£>2£) even more weakly with p 



lar to those obtained by Reinecke et al. 



Jl993| for the 
perpendicular diffusion, as shown in Fig. 3] as well. 

Given the uncertainties both in the actual magnetic 
turbulence evolution in the heliosphere with radia l dis- 



tance and latitude (see, e.g., Oughton et al. ( 20111) for a 



study in which they find a much more complicated ra- 
dial dependence of the slab variance) and their actual 
relation to perpendicular or even anisotropic perpendic- 
ular diffusion in connection with three-dimensional tur- 



bulence (Shalchi 2010 Shalchi et al. 2010) we stick, in 
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the following, with the empirical formulas of Eqs. (22) to 
( 24 ) for this principal study. 



4.3. The structure of the global diffusion tensor 

Employing the formalism described in Section 3 to cal- 
culate the global diffusion tensor k results in tensor el- 
ements Kij which are different from those 'traditionally' 
used, labeled nfj here, with i, j e {r, i?, ip}. The latter are 
derived follo wing the transformation presented by | Burger] 
et al. (2008), which for the Parker field is equivalent to 
the assumption that the local system can always be de- 
fined by t, ft — e§ x t, and b = e&, see Appendix A for the 
detailed transformation formulas. Both local systems are 
illustrated in Fig. 



The different behavior of the tensor elements Kij and 
nf, with latitude at a heliocentric distance of 5 AU and 
longitude ip — it /A is displayed in Fig. HI for the Parker 
and hybrid Fisk field, respectively. By definition, both 
formulations yield the same tensor elements in the eclip- 
tic plane, i.e. for ■& = ir/2, while for higher latitudes the 
differences become more and more pronounced. 

In the Parker case, e.g. the upper two rows in Fig. HJ 
the elements n rr , n rtp , and k w show roughly the same 
behavior for all latitudes. The strongest mixing of the 
local elements and k\ji occurs in k§§, so that the 
deviations for high latitudes are more pronounced. The 
main difference appears in the off-diagonal elements n r § 
and Ktfip, which are different from zero in the general 
case discussed here while they are equal to zero in the 
traditional approach. 

The differences between the hybrid Fisk field tensor el- 
ements (shown in the lower two rows in Fig. |4j) are simi- 
lar to those of the Parker field described above, although 
they show a more complicated $ dependence. Note that 
in the traditional formulation the off-diagonal elements 
K rl ? and Kdtp are already nonzero for the hybrid Fisk field 
and become larger in the new formulation. 

The choice of longitude is arbitrary for the Parker field, 
since it has no ip dependence. The Fisk field, however, 
has significant longitudinal variations, therefore, we show 
the ratios of the traditional and the new tensor elements 
for the hybrid Fisk field with longitude (Fig. [5| . It can be 
seen that for large heliocentric distances, the deviations 
between both formulations vary strongly, illustrated here 
for a heliocentric distance of r — 50 AU and a helio- 
graphic colatitude of d = ir/4. 

We emphasize again that in the case of isotropic per- 
pendicular diffusion (£ = 1), the traditional and the new 
formulations are identical for any given magnetic field 
with non-vanishing curvature. The differences between 
them scale with the perpendicular anisotropy £ (see Eq. 
24b. 



5. APPLICATION TO THE MODULATION OF COSMIC 
RAY SPECTRA 

To assess the impact of the new tensor formulation on 
CR modulation we employ a CR proton trans port model 
by solving the Parker equation (Parker 1965) 

dj_ 
dt 



v ■ (kvf) 



v/- 



(28) 



configuration space and p as momentum) . The solar wind 
velocity u s is radially pointing outwards with a constant 
speed of 400km/s and k is the diffusion tensor in the 
global frame for the Parker spiral magnetic field. This 
implies the Frenet-Serret trihedron of the form explicitly 
given in Appendix B. 

The solution is obtained via a numerical integration of 
an equiv alent system of stochastic differen tial equations 
(SDEs) flKopp et al.|[20TT| |Gardiner1[l994l ) 



dxi = Ai(xi) dt + ^^2 Bij(xi) dWj 



(29) 



for an ensemble of pseudo-particles (phase space ele- 
ments) with k = BB T and dW(t) = VdtN(t) where 

N(t) is a vector of normal distributed random numbers 
and Xi denotes the phase space coordinates. The stochas- 
tic motion dW(t) is often refereed to as Wiener process. 
The deterministic processes from Eq. ( 28 ) like the advec- 



tion with the solar wind flow and the adiabatic energy 

changes are contained in the generalized velocity A. We 
employ the time-backward Markov stochastic method, 
meaning that we trace back the pseudo-particles from a 
given phase space point of interest, until they hit the in- 
tegration boundary. The solution to the transport equa- 
tion ( |28[ ) is then constructed as a proper average over the 
pseudo-particle orbits. For details on the general method 
and the numerical sche me, especially in th e case of a gen- 
eral diff usion tensor, see Kop p et al. (2011 ), Strauss et al. 
( |2011a[) and |Strauss et al.| ( [21jTTby 

'1'he local interstellar spectrum (LIS) of protons jlis 
is assumed at a heliocentric distance of r = 100 AU 
as a spherically symmetric Dirichlet boundary condition. 
At the inner boundary of one solar radius r — Rq the 
pseudo-particles arc reflected, which is equivalent to a 
vanishing gradient in the CR density there. A standard 
representation of the proton LIS is given by 



3lis = Yl.lA[3{E kin +Q.bE ) 



-2.6 



to determine the differential CR intensity j(r,p,t) = 
p 2 f(r,p,t) (with f as the position in three-dimensional 



(30) 

and was taken from Reinecke et al. ( |1993 ). The proton 
rest energy Eq is equal to 0.938 and Ekin denotes the 
kinetic energy of a particle (both in units of GeV). 

The LIS and the resulting modulated spectra are 
shown in Fig|6]for both tensor formulations and for sev- 
eral heliocentric distances. The spectra for the new 
Frenet-Serret tensor are higher by up to 60% at low en- 
ergies for all heliocentric distances. This is due to the 
enhanced diffusive flux from the modulation boundary 
via an effective inward diffusion alo ng the polar axis. In 
the tensor formulation provided by Burger et al. (20081 
this diffusion (determined by n±2) cannot transport par- 
ticles from the boundary into the inner heliosphere, it 
merely distributes the particles on a shell of fixed he- 
liocentric distance. In the new tensor formulation exists 
thus a form of 'pseudo drift' produced by the off-diagonal 
tensor elements in the global frame, which were different 
or even equal to zero in the traditional formulation. This 
reduced modulation effect is relevant for higher energies 
at lower heliocentric distances, since the particles have 
more time to adiabatically cool (see the right panel of 
Fig.§. 
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Figure 4. The six independent tensor elements nf- for the "traditional" tensor formulation following Burger et al. (20081 (solid line) 



and the new «y using the Frenet-Serret trihedron (open circles) for a fixed radius of r = 5 AU, a longitude of ip = 7r/4 and for varying 
colatitude. The upper two rows show those for the Parker field while the lower two rows show those for the hybrid Fisk field. 
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hybrid Fisk field plotted against hcliographic longitude for a fixed 
heliocentric distance of r = 50 AU and a heliographic colatitude 
of i9 = 7r/4. The elements shown are k rr (red), k r $ (green), k rlp 
(blue), ktftf (violett), k$ v (brown), and k vv (black). 
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Figure 6. Modulated spectra for fully anisotropic diffusion of 
galactic protons for both tensor formulations. The left panel shows 
the resulting spectra for four heliospheric distances (1 AU, 25 AU, 
50 AU, 75 AU, from bottom to top) and the LIS modulation bound- 
ary at 100 AU (solid squares). The spectra, shifted in the plot by 
powers of 10 for clarity (note the resulting high energy offsets), 
converge to the LIS for high energies. While the symbols indicate 
the results from the new tensor formulation with the Frenet-Serret 
orientation, the lines are results from an analogous computation 
employing the 'traditional' two-angle transformation. In both cases 
is Kj_2 = 2/sxi- The right panel gives the relative deviations (nor- 
malized to the new results) of corresponding spectra from each 
other. The symbols are the same as in the left panel. 

6. CONCLUSIONS 

We have derived, in a global reference frame, the gen- 
eral form of the diffusion tensor of energetic particles in 
arbitrary magnetic fields. This new formulation partic- 
ularly includes the case of anisotropic perpendicular dif- 
fusion that arises from field line wandering or scattering 
due to turbulence and requires a determination of both 
principal (orthogonal) perpendicular diffusion directions. 
Unless the turbulence is non-axisymmetric, which ap- 
pears to be unlikely for the solar wind, the natural choice 
for these principal directions is the Frenet-Serret trihe- 
dron associated with the curvature and torsion of the 
magnetic field lines. 

After the derivation of the formulas for all tensor el- 
ements in dependence of the Frenet-Serret unit vectors, 



we have first quantitatively compared the results to those 
published previously for the example of the heliospheric 
magnetic field. For the latter we have discussed two well- 
established alternatives, namely the Parker field and the 
hybrid Fisk field. While the old and new tensor formu- 
lations coincide for the case of isotropic perpendicular 
diffusion, the more general case of anisotropic perpen- 
dicular diffusion cannot be treated consistently with the 
earlier approaches. This is manifest in significant differ- 
ences of corresponding tensor elements including addi- 
tional non-vanishing ones. 

Second, we have demonstrated the consequences of the 
new tensor formulation in application to the modulation 
of galactic cosmic ray proton spectra in the Parker he- 
liospheric magnetic field. Solving the cosmic ray trans- 
port equation with the method of stochastic differential 
equations allowed us to quantify the differences between 
the spectra resulting from both tensor formulations for 
the case of perpendicular diffusion with an anisotropy 
of £ = 2. We found those differences to amount up to 
60% at energies below a few hundred MeV. Given that 
we used for this first principal assessment an anisotropy 
that is moderate as compared to findings from detailed 
transport and modulation studies, the fluxes can be in- 
fluenced even more strongly and at even higher energies. 

Besides the fact that the above results indicate the ne- 
cessity to study the case of fully anisotropic diffusion in 
more detail within the framework of more sophisticated 
models of heliospheric cosmic ray modulation, they can 
furthermore be expected to be of importance for the par- 
ticle transport in complex galactic magnetic fields for 
which usually isotropic (scalar) diffusion has been con- 
sidered so far. 
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APPENDIX 

A. BURGER TRANSFORMATION FORMULAS 

The t ransformation formu las for the diffusion tensor 
given m |Burger et ah] ( |2008| ) read: 

,2 a i „„J//„ 2 ,t, i „ • 2 



ViS 
B 



"-00 — 

K B - 

K* = 



nj_ 2 sin C + cos C( K || cos * + K J-i sm 
sin C cos C(kii cos 2 * + sin 2 1 A7 — nj_ 2 ) 

— sin ^ cos 'J/ cos C( K \\ ~ K ±i) 

k± 2 cos 2 C + sin 2 cos 2 + k ±1 sin 2 *S>) 

— sin ^> cos "J sinC(ft|| — 



"ipxp 

with tan\f r 



km sin 2 + k ± i cos 2 * 



(Al) 



-By/ y jBt+ Bl and_tanC = B$/B r . 



Note that Kobylinski (2001) and Alania (2002) state a 
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different formula for "J, namely tan = —B v /B r . More- 
over, these formulas involve only two angles in contrast 
to the general case described with the matrix A in Eq. (1) 
in Section 2. As discussed in the text, these formulas in 
the given form can only hold for kh ^ k± 2 i n case of 
special magnetic fields with B$ = like that introduced 
by Parker. 

B. THE PARKER FRENET-SERRET TRIHEDRON 

Here, we derive the analytic expressions for the Frenet- 
Serret trihedron for the Parker case of the HMF. Reduc- 
ing Eqs. (18) to (20) to the Parker field by setting F s = 
we obtain 

eV-tanx^ ( B1 ) 



t = 



y/1 + tan 2 X 



-rsini?. The 



for the tangential vector, with tan^ = 

easiest way to derive the normal vector n is to calc ulate 
(t ■ V)f = kn (where k is the curvature, see Eq.(|12|) and 
to normalize appropriately. After some straightforward 
calculation, one arrives at 



_ £Ttan xe r + eL) + Fe# 

77, = — 

^/E 2 (ta,n 2 X + 1) +F 2 
where the abbrevations 



(B2) 



E = 



tan 2 x 



sin$ 



F = tan2xC » S ^ (B3) 
rsmv 



u s 1 + tan 2 x ' 

have been introduced. The binormal vector is now simply 
the cross product b = t X n, which yields 



b = 



-F(tan X e r + e v ) + 2?(tan 2 x + l)e# 
V / -F 2 (tan 2 X + 1) + # 2 (tan 2 x + I) 2 



(B4) 



Eqs.(Bl), (B2), and (B4) are the explicit formulas for 
the Frenet-Serret trihedron in the case of the heliospheric 
Parker field. Corres ponding but much longer expressions 
can, in principle, be obtained for the Fisk field as well. 
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